Reentrant behavior of divalent counterion mediated DNA-DNA electrostatic 

interaction 
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The problem of DNA-DNA interaction mediated by divalent counterions is studied using computer 
simulation. Although divalent counterions cannot condense free DNA molecules in solution, we show 
that if DNA configurational entropy is restricted, divalent counterions can cause DNA reentrant 
condensation similar to that caused by tri- or tetra-valent counterions. DNA-DNA interaction is 
strongly repulsive at small or large counterion concentration and is negligible or slightly attractive 
for a concentration in between. Implications of our results to experiments of DNA ejection from 
bacteriophages are discussed. The quantitative result serves to understand electrostatic effects in 
other experiments involving DNA and divalent counterions. 

PACS numbers: 87.19.xb, 87.14.gk, 87.16.A- 



The problem of DNA condensation has seen a strong 
revival of interest in recent years because of the need to 
develop effective ways of gene delivery for the growing 
field of genetic therapy. DNA viruses such as bacterio- 
phages provide excellent study candidates for this pur- 
pose. One can package genomic DNA into viruses, then 
deliver and release the molecule into targeted individual 
cells. Recently there is a large biophysic literature ded- 
icated to the problem of DNA condensation (packaging 
and ejection) inside bacteriophages dj. 

Because DNA is a strongly charged molecule in aque- 
ous solution, the process of ejection of DNA from bac- 
teriophages can be strongly influenced by the screening 
condition of the solution. By varying the salinity of so- 
lution, one can vary the amount of DNA ejected. Inter- 
estingly, monovalent counterions such as Na"*" have neg- 
ligible effect on the DNA ejection process In con- 
trast, multivalent counterions such as Mg+^, CoIIex+^, 
Spd^^ , or Spm"*"^ exert strong and non-monotonic effects 
[3|. There is an optimal counterion concentration, cz,q, 
where the least DNA genome is ejected from the phages. 
For counterion concentration, cz, higher or lower than 
czfi, more DNA is ejected from phages. The case of diva- 
lent counterions is more marginal. The non-monotonicity 
is observed for MgS04 salt but not for MgCl2 salt up to 
the concentration of lOOmM. 

The problem of DNA condensation by divalent coun- 
terions is a complex problem due to contributions from 
many physical factors. In the literature, most of the stud- 
ies dealing with this problem have focused on the ion- 
speciflc effects. For example, the hydration effects have 
been proposed to explain the above dependence on the 
type of divalent salts [3|. In this paper, we focus on role of 
non-specific electrostatic interactions between DNA and 
counterions. In a recent work we suggested that some 
aspects of DNA ejection in the presence of divalent coun- 
terions can be accounted for from the electrostatic point 
of view. Specially the strong, non-monotonic influence of 
divalent counterions on DNA ejection mentioned above 
is expected to have the same physical origin as the phe- 



nomenon of reentrant DNA condensation in free solu- 
tion [H, 0] . The fact that divalent counterions can have 
such strong influence on DNA ejection is not trivial. Un- 
like counterion of higher valences, Mg+^ counterions are 
known to not condense DNA 7|, or to condense them 
only partially in free solution [8|. However, DNA viruses 
provide a unique experimental setup. The constraint of 
the viral capsid strongly eliminates configurational en- 
tropic cost of packaging DNA. This allows divalent coun- 
terions to influence DNA condensation similar to that of 
tri- or tetra-valent counterions. In this paper, we use 
computer simulations to study the problem of DNA con- 
densation in the presence of divalent counterions. We 
show that indeed, if one includes only the non-speciflc 
electrostatic contribution, divalent counterions can in- 
duce DNA reentrant condensation like those observed for 
higher counterion valences. We offer an explanation for 
the discrepancy between DNA condensation in free solu- 
tion versus DNA condensation inside viruses. Our results 
show that, in addition to ion-specific effects, electrostat- 
ics exert a strong, non-negligible infiuence on qualitative 
and quantitative behaviors of this system. The results 
presented here can provide understanding of not only 
the electrostatics of DNA ejection problem, but also can 
serve as a starting point for investigating other systems 
involving DNA and divalent counterions where the phys- 
ical pictures are still not very well understood. 

We model the DNA bundle in hexagonal packing as a 
number of DNA molecules arranged in parallel along the 
Z-axis. In the horizontal plane, the DNA molecules form 
a two dimensional hexagonal lattice with lattice constant 
d (the DNA-DNA interaxial distance) (Fig. [T]). An 
individual DNA molecule is modeled as an impenetra- 
ble cylinder with negative charges glued onto it. The 
charges are positioned in accordance with the locations 
of nucleotide groups along the double-helix structure of 
a B— DNA. The hardcore cylinder has radius of 7A. The 
negative charges are hard spheres of radius 2A, charge 
— e, and lie at a distance of 9A from the DNA axis. This 
gives an averaged DNA diameter of Inm. The solvent 
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FIG. 1. (Color online) A DNA bundle is modelled as a hexag- 
onal lattice with lattice constant d. An Individual DNA 
molecule is modeled as a hard-core cylinder with negative 
charges glued onto it according to the positions of nucleotides 
of a B— DNA structure. 

water is treated as a dielectric medium with dielectric 
constant e = 78 and temperature T = 300°i^. The di- 
electric constant mismatch between water and DNA in- 
terior is neglected, and the cylinder only acts to prevent 
ion penetration. In our simulation, the positions of DNA 
molecules are fixed in space. This mimics the constraint 
on DNA configurational entropy inside viruses and other 
experiments of DNA condensation using divalent counte- 
rions. 

In the experiment of DNA ejection from bacterio- 
phages, there are both monovalent and divalent salts in 
solution. At very low concentration of divalent counteri- 
ons, DNA is screened mostly by monovalent counterions. 
To account for this limit, we include both salts in our sim- 
ulations. The mobile ions are modeled as hard spheres 
with unscreened Coulomb interaction (the primitive ion 
model). The radii of the colons and monovalent counte- 
rions are set to 2A. (For simplicity, we assume the two 
salts have the same colons.) The divalent counterions ra- 
dius is set to 2.5A. The interaction between two ions, i 
and j, with radii, Cij , and charges, qij, is given by 

U ^ f mj/enj if r,j >a,+ aj . 
\ 00 if Tij < ai + aj 

where = jr^ — r^j is the distance between the ions. 

The simulation is carried out using the periodic bound- 
ary condition. A periodic simulation cell with = 12 
DNA molecules in the horizontal (x, y) plane and 3 full 
helix periods in the z direction is used. The dimensions of 
the box are = 3d, Ly = 2V3d, and = 102A. The 
long-range electrostatic interactions between charges in 
neighboring cells are treated using the Ewald summation 
method. In Ref. [oj, it is shown that the macroscopic 
limit is reached when N > 7. Our simulation cell contains 
12 DNA helices, hence it has enough DNA molecules to 



eliminate the finite size effect. We did test runs with 1, 
4, 7, and 12 DNA molecules to verify that this is indeed 
the case. They are also used to check the correctness 
of our computer program by reproducing the results of 
DNA systems studied in Ref. |9|] in specific limits. 

In a practical situation, the DNA bundle is in equilib- 
rium with a water solution containing free mobile ions 
at a given concentration. Therefore we simulate the sys- 
tem using Grand Canonical Monte-Carlo (GCMC) sim- 
ulation. The number of ions are not constant during the 
simulation. Instead their chemical potentials are fixed. 
The chemical potentials are chosen in advance by simu- 
lating a DNA— free salt solution and adjusting them so 
that the solution has the correct ion concentrations. In a 
simulation, the ions are inserted into or removed from the 
system in groups to maintain the charge neutrality [lo| . 
Following Ref. , instead of using individual chemical 
potentials, /i+2, M-i-ii ^ind /^-i, for each ion species, we 
use only the combined chemical potentials, 

/^r" = M+2 + 2/i-i, M^i* =M+i+M-i, (2) 

in the Metropolis acceptance criteria of a particle inser- 
tion/deletion move. In this paper, we simulate DNA bun- 
dles at varying concentrations cz- Both /i^"/* and /x^^j* 
are adjusted so that the monovalent salt bulk concentra- 
tion, ci, in the DNA— free solution is always at 50mM 
(typical value of the DNA ejection experiment) and cz is 
at the desired value. Typical standard deviations in the 
final salt concentrations are about 10%. 

To study DNA— DNA interactions, we use the Ex- 
panded Ensemble method [§] to calculate the pressure 
of the DNA bundle. In this method, we calculate the 
difference of the system free energy at different volumes 
by sampling these volumes simultaneously in a simula- 
tion run. By calculating the free energy difference Afi 
for two nearly equal volumes, V and AV, we can cal- 
culate the total pressure of the system, P(T, V, {/Uj/}) — 
—Ail/AV (here {/ii/} is the set of chemical potentials of 
different ion species). The osmotic pressure of the DNA 
bundle is then obtained by subtracting the total pressure 
of the bulk DNA-free solution, Pt(r, V, {fJ-,^}), from the 
total pressure of the DNA system, Posm{T,V, {^i,}) = 
P{T,V,{fiA)-PbiT,V,{fi,}). 

In Fig. [5^, the osmotic pressure of the DNA bundle 
at different cz is plotted as a function of the interax- 
ial DNA distance, d. Because this osmotic pressure is 
directly related to the "effective" force between DNA 
molecules at that interaxial distance 0, Fig. also 
serves as a plot of DNA— DNA interaction. As one can 
see, when cz is greater than a value around 20mM, there 
is a short— range attraction between two DNA molecules 
as they approach each other. This is the well-known 
phenomenon of like-charge attraction between macroions 

[ll|. It is the result of the electrostatic correlations 
between counterions condensed on the surface of each 
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DNA molecule. The attraction appears when the dis- 
tance between these surfaces is on the order of the lateral 
separation between counterions (about 14A for divalent 
counterions) . The maximal attraction occurs at the dis- 
tance d ~ 28A in good agreement with various theoret- 
ical and experimental results 0, Il3 • For smaller d, the 
DNA-DNA interaction experiences a sharp increase due 
to the hardcore repulsion between the counterions. One 
also sees that the depth of attractive force between DNA 
molecules saturates at around —4 atm as cz increases. 
This saturation is easily understood. At small c^, there 
are both monovalent and divalent counterions present in 
the bundle. As Cz increases, divalent counterions replace 
monovalent ones in the bundle as the later ions are re- 
leased into the bulk solution to increase the overall en- 
tropy of the solution. However, charge neutrality condi- 
tion of the DNA macroscopic bundle and the hardcore 
repulsion between ions limit how many divalent counte- 
rions can be present inside the bundle. Once all mono- 
valent counterions are released into solution (replaced 
by divalent counterions), further increase in cz does not 
significantly change the number of divalent counterions 
in the bundle. This leads to the observed saturation of 
DNA— DNA short— range attraction with increasing cz- 

The strong influence of divalent counterions on DNA 
bundles can be seen by looking at the DNA-DNA "effec- 
tive" interaction at larger d. As evident from Fig. ^jji. for 
large d, at small cz DNA-DNA interaction is repulsive. 
As Cz increases, DNA-DNA interaction becomes less re- 
pulsive and reaches a minimum around lOOmM. As cz 
increases further, DNA-DNA repulsion starts to increase 
again. This non-monotonic dependence of DNA-DNA 
"effective" interaction on the counterion concentration 
is even more clear if one calculates the free energy of 
packaging DNA into bundles. This free energy is the dif- 
ference between the free energy of a DNA molecule in a 
bundle and that of an individual DNA molecule in the 
bulk solution. Per DNA nucleotide base, this free energy 
is given by: 

MDNA(^) = (i/LzN) f Posm{d')dV{d'), (3) 

J oo 

here I = 1.7 A is the distance between DNA nucleotides 
along the axis of the DNA, and V{d) — L^LyL^ — 
6^/3Lzd^ is the volume of our simulation box. The re- 
sult for optimal bundle lattice constant 
d — 28A is plotted in Fig. [2b as function of the cz 
[ist . Once again, there is an optimal concentration, 
czfi — lOOmM, where the free energy cost of packaging 
DNA is lowest. It is even negative indicating the ten- 
dency of the divalent counterions to condense the DNA. 
At smaller or larger concentrations of the counterions, 
the free energy cost of DNA packaging is higher. 

This reentrant behavior of DNA interaction can be un- 
derstood. At large separations, the distribution of coun- 
terions in the bundle can be considered to be composed 
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FIG. 2. (Color online) a) The osmotic pressure of the DNA 
bundle as a function of the interaxial DNA distance d for 
different divalent counterion concentration cz shown in the 
inset, b) The numerically calculated free energy of packaging 
DNA molecules into hexagonal bundles as a function of the 
divalent counterion concentrations. 



of two populations: condensed layers of counterions near 
the surfaces of the DNA molecules and diffuse layers of 
counterions further away. It is reasonable to expect the 
thickness of the condensed counterion layer to be on the 
order of the average lateral distance between counterions 
on the DNA surface (w 14A) . So for d > 34A, both coun- 
terion populations are present and one expects DNA- 
DNA interaction to be the standard screened Coulomb 
interaction between two charged cylinders with charge 
density rj* . The qualitative dependence of 77* on cz can 
be obtained by plotting the local colon concentration c_i 
as a function of distance from DNA axis (Fig. [3]). At 
low Cz, c_i decreases as d decreases from 00 suggesting 
1]* is negative (undercharged DNA). At high cz, c_i in- 
creases as d decreases until the condensed counterion lay- 
ers start to overlap at d « 34A. This shows that 77* is pos- 
itive (overcharged DNA). In both cases, DNA repulsion 
is strong. For an intermediate value of cz, ?7* — 0, DNA 
is almost neutral, and the repulsion is weakest. Further- 
more, the like-charge attraction among DNA molecules 
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FIG. 3. (Color online) Local concentration of colons as a 
function of distance from the axis of a DNA in the bundle for 
different divalent counterion concentrations, for d = 50A. 



mediated by the counterions is dominant in this con- 
centration range, causing the electrostatic packaging free 
energy to become negative. 

Figure gives a value of — O.OOl/csT/base for the 
short— range attraction among DNA molecules at the op- 
timal concentration. This is slightly less negative than 
revious theoretical fit of viral DNA ejection experiments 
We believe this small difference is due our choice of 
the system's physical parameters such as ion sizes [14 1. 
The azimuthal orientation correlations of DNA [l^ are 
another omission in our study. Relaxation along this de- 
gree of freedom can further lowering energy of the sys- 
tem. The non-electrostatic (such as van der Waals) in- 
teractions at small d can also enhance DNA attraction. 
On the other hand, dielectric constant mismatch between 
water and DNA interior could push the condensed ions 
away from DNA interior and lower the attraction energy. 
However, these effects are minor at large DNA— DNA 
separations, thus do not change the qualitative reentrant 
condensation picture. More comprehensive studies that 
take these effects into account are the subjects of our fu- 
ture works. Nevertheless, the value range of — O.OOI/cbT 
obtained in this paper is significant. It explains why diva- 
lent counterions exert strong effect on DNA ejection from 
virus but are not able to condense DNA in free solution. 
This value corresponds to an attraction of a fraction of 
—ksT per one persistence length (~ 300 bases). This is 
too small to overcome thermal fluctuation of DNA (about 
one ksT per persistence length), thus cannot condense 
them. Only inside the confinement of the viral capsid, 
where DNA configuration entropy is strongly suppressed, 
can divalent counterions cause strong influence. The non- 
monotonic behavior described above has the same physics 
as the phenomenon of reentrant DNA condensation by 
counterions 0, [3l of high valences. In this paper, we 
demonstrate clearly that it can happen to divalent coun- 
terions if DNA configuration entropy is restricted. This 



correlates well with experimental data of DNA ejection 
from bacteriophages. We should mention here that DNA 
condensation by divalent counterions has also been ob- 
served in another environment where DNA configuration 
is constrained, namel y th e condensation of DNA in two 
dimensional systems [17[. This fact once again strongly 
supports our argument. 

In conclusion, in this paper, we use a computer sim- 
ulation to study the electrostatics of DNA condensation 
in the presence of divalent counterions. The entropy of 
DNA configure fluctuation is suppressed in simulation. 
Such study can be applied directly to the experimen- 
tal problem of DNA ejection from bacteriophages where 
DNA condensed in a strongly confined environment. Our 
results show that, even at the level of non-specific elec- 
trostatic interaction, divalent counterions can strongly 
influence DNA ejection. This potentially opens up an ad- 
ditional degree of freedom in controlling bacteriophages 
for the general purpose of gene therapy or viral diseases 
treatment. Beyond the scope of DNA ejection experi- 
ments, we believe the quantitative results of our paper 
can be used to understand many other experiments in- 
volving DNA and divalent counterions. 
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